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We use the distribution, and particularly the skewness, of high redshift type la supernovae bright- 
nesses relative to the low redshift sample to constrain the density of macroscopic compact objects 
(MCOs) in the universe. The data favors dark matter made of microscopic particles (such as the 
LSP) at 89% confidence. Future data will greatly improve this limit. This constraint is valid for 
a range in MCO mass from 10 -2 Mq to 10 10 Mq. Combined with other constraints, MCOs larger 



than one tenth the mass of Earth (^ 
matter. 
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can be eliminated as the sole constituent of dark 



The recent observational success of the cold dark mat- 
ter (CDM) model puts only loose constraints on the 
properties of dark matter itself. The requirements for 
a cold dark matter candidate can be satisfied by a mi- 
croscopic, weakly interacting particle or by macroscopic 
objects that do not produce significant amounts of ob- 
servable radiation through astrophysical processes, such 
as primordial black holes (PBHs). Direct observational 
constraints are needed. 

There have been a number of previous constraints put 
on the density of massive compact objects (MCOs) from 
a variety of observations. The rate of gravitational mi- 
crolensing of stars in the Small and Large Magellanic 
Clouds puts a limit on the mass fraction of MCOs in 
our galaxy's dark matter halo of < 40% with masses 
between 10 -7 and 10 M Q assuming a "standard" halo 
model P, 0, Microlensing toward M31 has thus far 
been inconclusive 0, [B[ . From the effect of gravitational 
lensing on the equivalent widths of quasar broad lines Q 
placed a limit of fl mco < 1 for 10~ 3 M < M < 300 Mo- 
using Very Large Baseline Interferometry (VLBI), pj 
looked for radio QSOs that are gravitationally lensed into 
multiple images. They put an upper limit of f2 mco < 0.01 
for very massive (10 6 - 10 8 M ) MCOs. There are also 
arguments from stellar dynamics for why such high mass 
MCOs cannot make up a significant fraction of the Milky 
Way's dark matter halo @, Q . Although these arguments 
eliminate the possibility of all the dark matter being in 
this form it is difficult to be more quantitative since we 
do not have detailed information on the initial population 
of globular clusters or the velocity dispersion of the disk 
stars. All of these constraints leave a large open window 
between ~ 100 M and 10 5 M which this study seeks 
to eliminate. 

Type la supernovae are found empirically to be good 
standard candles with luminosities large enough to be ob- 
servable at cosmological distances. They have been used 
to measure the expansion rate of the universe as a func- 
tion of time and through this the densities of matter and 



dark energy 1(| 11, and references therein]. The distri- 
bution of supernova brightnesses will also be systemati- 
cally affected by the gravitational fields along the line of 
sight. In the most extreme case where all the matter in 
the universe is in macroscopic objects, the magnification 
distribution will be highly skewed with most high red- 
shift supernovae appearing less bright than they would 
in a perfectly homogeneous universe and a small number 
being highly magnified [H, 13, [3, EH- In this paper 
we will use this effect to measure what fraction of dark 
matter can be made of MCOs. The systematic skew in 
the distribution of luminosity distances can also affect 
the implied best-fit cosmological parameters. The cos- 
mological constant is likely to be overestimated if MCOs 
exist in large numbers but are not included in the data 
analysis. We show how significant this mis-measurement 
could be. 

Regardless of what the basic unit of dark matter is, 
SNe will be weakly lensed by the clumping of the dark 
matter and baryons into galaxies, dark matter halos and 
larger scale structures [3, H3] • However this effect is 
relatively weak and as yet undetected [ll| (see [l9| for a 
contrary opinion). We choose the conservative option of 
not including this type of weak lensing in our analysis. 
We derive an upper bound on the density of MCOs and 
including clumping would only reduce this bound at the 
expense of further assumptions. 

With the use of the Sachs optical scalar equations [2(| 
the cross section of a beam of light can be calculated 
as a function of the affine parameter along a fiducial 
light path, A. The beam is affected by both the local 
mass-energy density that the beam passes through, rep- 
resented by the Ricci tensor, and by nonlocal, inhomo- 
geneous shear contributions coming from mass outside of 
the beam which are proportional the the Weyl tensor. Ig- 
noring the Weyl contribution, the angular size distance 
can be define as the solution to the following equation 
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It has been assumed that the light beam passes through 
a fluid of density p(A) and negligible pressure (p << 
pc 2 ). The affrne parameter is related to the red- 
sh ift by dX = cH^jl + z)~ 2 E(z)- 1 ^ whe re E(z) = 

y/n m (i + z y + n A + (i - n m - n A )(i + z) 2 . 

The standard angular size distance calculated with the 
Friedman-Robertson- Walker (FRW) metric (the angular 
size distance in a perfectly homogeneous universe), here 
denoted by D(z, f2 m , ^a), is a solution to ([T]) in the spe- 
cial case of p(z) = p(z) where p(z) is the average density 
of the universe. The angular size distance appropriate for 
a light beam that travels through empty space or through 
a density of matter that is less than the average density 
will be larger than the FRW distance. The angular size 
distance between two redshifts z\ and z^ for a beam that 
passes through a fixed fraction of the average density, 
/ mco , will be denoted D(z\, z 2 > fmco, ^, ^a)- We call 
this the unfilled beam distance when f mco < 1. 

If a point-like source is being lensed by a single point- 
like mass, the source position on the lens plane, y, 
can be expressed in terms of the magnification y = 
V2R e (z s J mco ,n m ,{l A )^p(p 2 - lp72 _ l gj. There 
are two images of the source and p is the sum of their 
magnifications. The Einstein radius is given by 



Re(z , Z S: /mco > ^mj ^a) 
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The distances that appear in the Einstein radius are the 
unfilled beam distances and not D(z). The magnifica- 
tion probability distribution function (pdf) is easily cal- 
culated 



P{p) dp = 2tt 
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dz M dp 
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where \ is the coordinate distance in the radial direction. 
The normalization, r (z s , / mco , J7 m , J7a), is often called 
the optical depth. Note that the optical depth depends 
only on the average mass density of MCOs and not on 
the mass of the individual MCOs. It is also independent 
of the Hubble parameter since D cx Hg 1 . 

The magnification pdf for one lens is not sufficient for 
an ensemble of point-like lenses. There is no way of ana- 
lytically calculating the full pdf for an ensemble of lenses. 
However, there are several requirements that the full pdf 
must fulfill. One is that the average angular size dis- 
tance to objects at a fixed redshift be the same as it 



would be in a homogeneous universe [23j. This is the 
requirement that (p)D(z)~ 2 = D(z)~ 2 . We also expect 
that for p < 1 the pdf is zero or at least very small. 
This is because when the optical depth is small (as it 
is for all the cases considered here) the p ~ 1 cases oc- 
cur for lines of sight that pass well away from any lens 
where there will be only one significant image. The weak 
lensing limit (/ drD(0, z)D(z, z s ) VV<K r ) < c 2 D(0,z s ), 
where 4» is the Newtonian potential and the integral is in 
the radial direction) will apply in these cases and in this 
limit there must be at least one image with p > 1 [24| . 

In addition, we expect that when the optical depth is 
small the true pdf will have the same form at large p as 
([3]) because in these rare cases the line of sight passes 
very close to one lens and is unlikely to pass very close to 
two or more. Numerical simulations give a magnification 
pdf with a high p tail very similar to ([3]), as expected, 



and a strong, but not singular mode at p = 1 12, 3, 14 1. 
In this paper, to estimate the mag nification pdf we use 
an analytic form advocated by [12j: 



P(p; S, T eS )dp 



2t, 



off 



l_ oxp (i-^) 




3/2 



dp , p > 1 
, p< 1 
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We fix the two free parameters, 5 and T e g, by applying 
the average magnification constraint and the normaliza- 
tion constraint numerically. We have experimented with 
other forms for the pdf and find that any pdf that satis- 
fies the expectations given above and is monotonic and 
reasonably smooth will give very similar results for the 
density of MCOs. 

We use two sets of p ublicly available SN data. Firstly, 
the Riess et al. [llj (here after RST) "gold" sample 
containing 157 type la SNe. There is an empirical re- 
lationship between the width of the SN's light curve and 
the peak luminosity of the SN. They derive template 
light curves in different passbands from an independent 
sample of SNe [25| and then fit the light curves of the 
high redshift SNe to these templates to determine a cor- 
rected peak luminosity for each SN. Because the correc- 
tion method is derived from a low redshift sample of SNe 
that should not be affected by gravitational lensing it can 
be justifiably used to look for lensing in the high redshift 
sample. 

Secondly, the Supernova Legacy Survey (SNLS) [26[ 
uses 44 "nearby" SNe (0.015 < z < 0.125) from the litera- 
ture, some of which are included in the RST gold sample, 
plus 73 of their own high redshift SNe (0.249 < z < 1.01). 
SNLS makes the peak-luminosity correction by the sim- 
ple relation m \, = m* B + a(s — 1) — (3c where s is the 
"stretch", a measure of the light curve width, c is the 
color (= B — V + 0.057) and m* B is the observed peak 
magnitude in rest frame i?-band. These three parameters 
are fit to each observed light curve. The parameters a, 
(3 and the absolute peak magnitude of the SNe are fit to 
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the whole data set along with the cosmological model. In 
what follows we refit these parameters since the fits done 
by SNLS were done without including the possibility of 
gravitational lensing. 

When gravitational lensing is included, the probabil- 
ity distribution for the magnitude of a SN is not Gaus- 
sian and x 2 is n °t a useful statistic. We instead adopt a 
Bayesian approach to constrain the parameters. For each 
SN we use a likelihood function of the form 



I-oo^f* P H (P\ z ») P m{™>ob, Cob, A*) 



(5) 



where m b is the observed, light-curve-corrected magni- 
tude, G b is the reported observational error. 

It is probable that the intrinsic distribution of cor- 
rected peak luminosities is not Gaussian distributed. If 
the intrinsic distribution has a skew it could be mistaken 
for a lensing effect (if negative) or it could obscure a 
lensing signal (if positive). To account for this we allow 
the intrinsic distribution of absolute magnitudes to have 
a redshift independent skewness. This is implemented 
through modeling P rn with a truncated Edgeworth ex- 
pansion 



Pn 



-x 2 /2 
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(m o6 - 5 log(£> L ) + 2.5 log(^) - M a ) 
x —— — 

V°l + °Tb 

where K3 is the skewness. The lower redshift SNe, which 
are not lensed, constrain K3. This gives a total of 6 free 
parameters for the RST SNe. For SNLS SNe there are 
two additional parameters (a and 0). 

Recent CMB observations and galaxy surveys strongly 
indicate that the universe is flat (O m + Q\ = I) or very 



close to it [27[ ■ In this paper we will restrict ourselves to 



the flat case. We also use an independent constraint on 
f2 m from WMAP and the HST Key project measurement 
of the Hubble parameter as reported in [28| . Specifically 
Q m = 0.28 ± 0.04. This constraint is implemented as a 
Gaussian prior in the likelihood function. For all other 
parameters we use uniform priors. 

To find the confidence regions for all the parameters 
we perform a Monte Carlo Markov Chain (MCMC) cal- 
culation. This allows use to calculate a constraint on 
/mco (or any other parameter) that is marginalizing over 
all the other parameters. In each MCMC at least 60,000 
combinations of the parameters are tested. Two sets of 
calculations are done - one in which / su b is allowed to 
take any value between and 1 and a second set in which 
/sub is either or I. 

The recovered values for a, (3, the a m 's and the M a 's 
are all in good agreement with what was reported in the 
data papers. Considering first the comparison between 
the two hypotheses / su b = and 1, we find that the SNLS 
data proffers / su b = at 89% confidence. The RST data 



1 .0 
0.8 
I0.6 

E 

2. 0.4 

0.2 

0.0 
10'° 

10 s 

a 10 6 

D 

E 10" 

8 10 2 

5 10° 

10" 2 
10" 4 





5% chance of significant time delay 



5% chance of time dependence 




.0 0.2 0.4 0.6 0.1 



1 .0 



FIG. 1: On top is the Bayesian probability that the MCO fraction 
is less than / mco . The top curve is for the SNLS SN sample and 
the bottom is for the gold sample of [ill ]. The dotted line is F(< 
/mco ) = /mco ■ In the bottom panel the two solid curves show where 
the approximation of point-like source and zero time delay begin to 
breakdown and time-dependent lensing effects are expected. The 
hashed region is ruled out at 90% c.l. by the SN data plus the 
WMAP constraint Oil Qrn • 



proffers / su b = 1 at 84% confidence. These seem in con- 
tradiction until the K3 distributions are considered. For 
the RST data there is a strong positive correlation be- 
tween K3 and / su b while for the SNLS data they are not 
correlated. If K3 is fixed at the RST data shows no 
preference between / su b = and 1 while the preference 
persists in the SNLS data. A maximum likelihood com- 
parison test shows that including K3 ^ does not im- 
prove the fit significantly enough that it is required by 
the data. For this reason we do not think that the RST 
data contains enough information on the skewness of the 
brightness distribution to be useful for constraining / su t>- 
This difference is probably a result of the different meth- 
ods for correcting the peak luminosity. 

Constraints on continuous values for / mco are shown 
in the top panel of figure [TJ The SNLS puts a Bayesian 
upper limit on / mco of 88% at a 95% confidence level. 

To see if including the possibility of MCOs in the anal- 
ysis changes the best fit VL m we repeat the calculations 
without the prior constraint on fl m . The best fit £l m 
for the SNLS data does not change significantly. There 
is a larger although still small effect for the RST data 
(n m = 0.29 ± 0.03 to 0.31 ± 0.05 when the prior is re- 
moved) . 

So far the SNe have been treated as perfect point 
sources. Since a SN's photosphere does have a finite size 
there is a lower mass bound where this approximation 
fails to hold. This comes about when the source size, 
projected onto the lens plane l s , is of order R e or smaller 
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and the impact parameter y is of order the sources size 
or smaller 29|, [3(|. We calculate the magnification nu- 
merically approximating the surface brightness profile as 
a uniform disk (the photosphere is expected to be opti- 
cally thick during these stages of the explosion) and take 
a deviation of 0.5 mag from the point source value to 
be significant. This requirement puts a upper limit on y 
as a function of l s /R e for a lens that could violate the 
point source approximation. Using this we calculate the 
probability that there is a lens that meets these condi- 
tions in the same manner that t(z s , / mco ) was calculated 
previously. 

The expansion rate of type la SN photospheres is 
Vphot = 1.0 - 1.4 x 10 4 km s" 1 31]. We take the typ- 
ical period in which a SN is observed to be 20 days. If 
the SN reaches a size where the point source approxima- 
tion is violated during this period the SN's light curve 
will be significantly affected. In the bottom panel of fig- 
ure Q] the curve is shown where there is a 5% chance that 
more than 1% of the SNe are affected in this way. The 
SN redshift distribution is taken to be the same as in the 
RST sample. This can be considered a conservative lower 
mass limit on where our constraints on / mco are valid. 

The magnification can also become time-dependent 
when the MCO mass scale is very large because of the 
time delay between the gravitationally lensed images. If 
this time delay is large enough compared to the width 
of the peak in the light curve and the images have sim- 
ilar magnifications then the light curve will be affected. 
The rise time for a type la SN is much shorter than the 
fall time so the fractional change in the brightness, 8f, 
reaches a maximum when the trailing image reaches the 
maximum. For the purposes of calculation we approx- 
imate the light curve as exponentially falling from its 
maximum with a time scale of 20 days. We consider a 
significant change in the light curve to be Sf = 0.06 cor- 
responding to 3 magnitudes below the peak magnitude. 

For a lens to cause a change in the light curve larger 
than Sf the SN must be within a range in x where the 
outer limit is set by when the magnification of the trailing 
image is too small to cause a significant change and the 
inner radius is set by the requirement that the time delay 
be large enough relative to the fall time of the SN. For 
a small enough mass it is not possible to satisfy both of 
these requirements at the same time. The probability of 
a lens being in this range can be calculated in a similar 
way to the calculation of t(z s , f mco ). The curve where 
there would be a 5% chance of having a significant change 
in the light curve of at least one of the SNe in the RST 
sample is shown in the bottom panel of figure [1] In effect 
this probability is independent of / mco except at very 
small values of / mco - This curve marks a formal upper 



mass boundary of where our / mco limit is valid. 

If M is outside these upper and lower bounds some 
of the SN light curves should show significant deviations 
from the standard shape. Unfortunately there are not 
good limits on time dependent deviations in the data 
because of irregular sampling and intrinsic variations in 
the light curve shapes. 

It has been shown that, at 95% confidence, no more 
than 88% of the dark matter in the universe can be 
in compact massive objects with masses greater than 
10 -2 M , closing a gap in the allowed mass for dark 
matter particles. Future data will undoubtedly improve 
these constraints. Gravitational lensing cannot signifi- 
cantly change the conclusions drawn from the high red- 
shift type la supernova surveys, but it will be important 
to account for it in future high precision measurements 
of the cosmological parameters. 
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